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ABSTRACT 


The  fracture  dynamic  and  crack  arrest  responses  of  modified  compact 
tension  specimen  (M-CT)  machined  from  Homalite-100  and  polycarbonate  sheets 
were  studied  by  dynamic  photoelasticity,  dynamic  finite  element  analysis  and 
streaking  photography.  In  contrast  to  the  results  of  a previous  study  involving 
a mild  steel  M-CT  specimen,  substantial  dyn2unic  effects  were  observed  during 
crack  propagation  in  the  Homalite-100  and  polycarbonate  M-CT  specimens.  Although 
the  crack  arrest  toughnesses,  Kjg,  were  within  10  percent  of  the  corresponding 
static  stress  intensity  factor  at  crack  arrest,  their  values  were  about  80  percent 
and  50  percent  of  the  corresponding  fracture  toughness,  Kj^.,  of  Homalite-100  and 

polycarbonate,  respectively. 

INTRODUCTION 

In  a recent  paper,  two  of  the  authors  used  a dynamic  finite  element  code  to 
compute  the  dynamic  fracture  toughness  of  a fracturing  transverse  wedge-loaded, 
modified  compact  tension  (M-CT)  specimen  machined  from  AISI  1018  steel  [1].  The 
dynamic  finite  element  code,  HONDO  [2],  was  used  in  its  "generation  phase"  [3] 
where  the  crack  was  driven  by  the  experimentally  determined  crack  motion  and  the 
associated  dynamic  fracture  toughness,  Kjp,  was  calculated.  This  singular  result 
involving  a very  ductile  material  with  a notched  brittle  weld  starter  crack  [4], 
lead  to  the  conclusion  that  little  difference  between  the  dynamic  and  static 
stress  intensity  factors  existed  in  the  particular  M-CT  specimen  analyzed.  This 
result  was  not  only  in  disagreement  with  all  previously  obtained  results  for 
different  specimen  geometries  [5]  but  also  contradicted  the  analytical -experimental 

results  obtained  by  Hahn  et  •alv  {69 'a' similar  M-CT  specimens  but  with  differ- 
ent  crack  starter.  Althou^.a  va4id  strain  fracture  toughness,  Kj^,  for 

the  mild  steel  M-CT  specimen  analyzed  in  Reference  [2]  was  not  available,  the 
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calculated  crack  arrest  stress  intensity  factor,  Kj^,  of  approximaiely 
88  MPg/lir  (80  ksi  /in  ) would  have  been  considerably  lower  than  such  fracture 
toughness  and  was  not  consistent  with  the  photoelastic  results  obtained  for 
Homalite-100  M-CT  specimens  by  T.  Kobayashi  et  al . [7]  where  Kj^  was  nearly 

equal  to  the  Kjc  = 445  KPa/lT  (405  psi  /Tn  ) . 

Since  similar  differences  between  Kjq  and  Kj^  were  observed  in  dynamic 
tear  (DT)  specimens  machined  from  the  relatively  brittle  Homalite-100  and  the 
ductile  polycarbonate  plates  [8,9],  it  is  thought  that  the  influence  of  duct- 
ibility  could  be  delineated  if  comparative  fracture  dynamic  studies  were 
conducted  on  M-CT  specimens  machined  from  Homalite-100  and  polycarbonate  plates. 

As  a result,  a combined  experimental -numerical  analyses  of  M-CT  specimens  mach- 
ined from  Homalite-100  and  polycarbonate  were  conducted  and  reported  in  this 
paper. 

The  M-CT  specimens,  as  shown  in  the  legend  of  Figure  1,  are  full-size  models 
of  the  dynamic  fracture  specimen  being  investigated  in  a current  ASTM  E24.03.04 
Subcommittee  on  Dynamic  Testing,  Dynamic  Initiation-Crack  Arrest  Task  Group  [10]. 
The  experimental  and  numerical  procedures  used  in  this  study  were  the  now  popular 
dynamic  photoelasticity  [7,8,9]  and  the  dynamic  finite  element  method  [8,9]  used 
in  its  generation  phase  [3],  respectively.  In  three  experiments,  crack  velocity 
measurements,  which  were  in  the  past  obtained  by  discrete  crack  length  measure- 
ments and  the  timing  marks  from  a Lite-Mike,  were  also  obtained  from  continous 
crack  length  recording  using  a streaking  camera. 


EXPERIMENTAL  PROCEDURES 


Dynamic  Photoelasticity 


The  16-spark  gap  Cranz-Schardin  camera  and  the  associated  dynamic  photo- 
elasticity system,  which  was  originally  developed  by  Riley  and  Dally  [11],  has 


been  discussed  In  many  previous  publications  and  thus  will  not  be  repeated  here. 
Figure  2 shows  two  typical  dynamic  photoelastic  patterns  surrounding  a running 
crack  tip  In  a polycarbonate  M-CT  specimen.  The  two  polycarbonate  M-CT  specimens 
of  6.4  mm  (1/4  In.)  thickness,  which  were  analyzed  by  dynamic  photoelasticity, 
were  annealed  at  160*’C  overnight  to  eliminate  residual  stresses.  A starter 
crack  of  approximately  2.5  mm  (1/4  In.)  length  was  sawed  and  chiseled  from  the 
tip  of  the  machined  notch.  The  blunt  starter  crack  Initiated  crack  propagation 
at  a relatively  high  crack  Initiation  fracture  toughness,  Kjq,  and  thus  propagated 
the  crack  nearly  through  the  entire  width  of  the  specimen.  The  average  mechanical 
and  optical  properties  used  In  photoelastic  data  reduction  as  well  as  In  dynamic 
finite  element  analyses  are  identical  to  those  listed  in  Reference  [9]. 

Errors  in  Kjp  Determination 

When  dynamic  photoelasticity  is  used  for  fracture  dynamic  analysis,  the 
transient  dynamic  isochromatic  lobes  surrounding  the  running  crack  tip  must  be 
related  to  an  instantaneous  dynamic  stress  intensity  factor.  Following  the 
original  suggestion  by  Irwin  [12],  normally  this  conversion  is  made  by  using 
either  one,  two  [13,14]  or  multiple  terms  [15,16]  in  the  static  crack-tip 
stresses  of  Williams  eigenfunction  [17]. 

Errors  involved  in  using  the  above  static  near-field  state  of  stress  were 
later  assessed  by  Kobayashi  and  Mall  [18]  who  used  the  dynamic  counterpart  [19] 
of  Williams  stress  function  to  show  that  overestimations  of  10  percent  or  more 
In  dynamic  stress  Intensity  factor  were  possible  at  a relative  slow  crack  velocity, 
a,  of  15  percent  of  the  dllatatlonal  wave  velcoty,  I.e.  0.15c^,  and  that  such  error 
Increased  with  the  use  of  larger  dynamic  Isochromatic  lobes.  Reference  [18]  also 

showed  that  stress  waves  propagating  In  the  vicinity  of  the  crack  tip  could 
distort  the  dynamic  Isochromatic  lobes  and  could  i.iduce  additional  errors  In 
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dynamic  stress  intensity  factor  determination.  Such  stress  waves  could  be  a 
visible  rectangular  pulse  as  recorded  photoelastically  by  Wade  and  Kobayashi  [20] 
or  an  innocuous  ramp  pulse  which  is  about  to  impact  the  propagating  crack  tip. 

As  a result  of  possible  compounded  errors  involved  in  the  use  of  statis 
stress  field  to  characterize  a dynamic  phenomena  and  from  the  not-so-apparent 
stress  wave  interactions  with  the  propagating  crack  tip,  the  authors  have  used 
the  smallest  visible  isochromatic  lobe,  prefereably  within  2.5  mm  (0.1  in.) 
distance  from  the  moving  crack  tip,  to  extract  the  dynamic  stress  intensity 
factor  at  higher  crack  velocities.  Such  size  restriction  on  the  permissible 
isochromatic  lobe  unfortunately  taxes  the  exprimental  accuracies  in  determining 
the  size  of  the  isochromatic  lobe  as  well  as  the  instantaneous  location  of  the 
crack  tip.  This  limited  resolution  in  Kjq  determination  via  dynamic  photo- 
elasticity, which  is  estimated  to  be  at  the  best  +5  percent,  is  akin  to  the 
corresponding  limitation  in  dynamic  finite  element  analysis  when  used  for 
dynamic  fracture  analysis. 

Crack  Velocity  Measurements 

In  a previous  paper  [18],  the  authors  discussed  the  experimental  errors 
involved  in  measuring  the  crack' tip  motion  and  the  need  to  smooth  the  raw  data 
to  confine  the  experimental  scatter  in  the  Kjq  versus  a relation  of  Reference 
[21].  Recent  numerical  experimentation  [22]  using  an  upgraded  fracture  mechanics 
subroutine  in  a dynamic  finite  element  code  shewed  that  slight  perturbations  in 
the  crack  tip  motion,  which  resulted  in  mild  oscillations  in  the  crack  tip 
velocities,  could  generate  significant  osciallations  in  the  calculated  Kjq. 

In  order  to  determine  the  existence  or  lack  of  existence  of  crack  velocity 
variation  during  dynamic  crack  propagation,  crack  velocities  were  measured  in 
two  Homalite-100  (thickness  9.5mm)  and  one  polycarbonate  M-CT  specimens 
(thickness  6.4mm)  using  a Beckman  Whitely  Model  318  streaking  camera. 
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Figure  3 shows  schematically  the  experimental  setup  as  well  as  a typical  streak- 
ing photograph  of  a fracturing  polycarbonate  M-CT  specimen.  Simultaneous  dynamic 
photoelastic  recording  was  not  possible  because  the  stray  light  from  the  light 
source  for  the  streaking  camera  Interfered  with  the  Crane-Shardin  camera  system. 
The  apparent  high  Initial  crack  velocity  from  the  streaking  photograph  of  Figure 
3 Is  due  to  the  lack  of  an  adequate  pre-triggering  system  for  the  light  source. 

I The  estimated  crack  velocity  at  the  onset  of  rapid  crack* propagation  was  thus 

I extrapolated  from  the  steady  state  crack  velocity  as  marked  In  the  streaking 

photography  of  Figure  3. 

Figures  4 and  5 show  typical  crack  velocity  relations  generated  from  the 
crack  position  versus  time  relations  obtained  by  the  streaking  photographs. 

The  crack  velocities  in  the  two  Homallte-100  and  one  polycarbonate  M-CT  specimens 
exhibited  little  change  during  much  of  the  crack  propagation. 

DYNAMIC  FINITE  ELEMENT  ANALYSIS 

The  dynamic  finite  element  code  which  was  Initially  [2]  introduced  for 
dynamic  fracture  analysis  in  the  "generation  mode"  has  undergone  substantial 
• changes  in  the  past  four  years.  An  Improved  crack  tip  release  mechanism  for 
rapid  crack  propagation  has  been  developed  and  an  updated  plane  stress  algorithm 
for  computing  dissipation  energy  at  the  crack  tip  has  been  Incorporated.  After 
many  numerical  experimentations,  a linearly  varying  crack  tip  nodal  release 
force  was  found  to  adequately  simulate  a more  gradual  transition  of  the  crack 
tip  movement  to  Its  adjacent  node  [23].  An  additional  Improvement  made  for  this 
study  Is  an  Iteration  algorithm  during  each  built-in  time  Increment  of  HONDO  [2] 
to  match  the  applied  nodal  force  with  the  nodal  force  calculated  from  the 
Incremental  change  In  nodal  velocity  In  this  explicit  dynamic  finite  element 
code.  Typically,  satisfactory  convergence  of  this  Iteration  scheme,  as  shown 
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In  Figure  6,  is  obtained  on  the  average  within  three  iterations  and  thus  the 
computational  efficiency  of  HONDO  is  still  preserved  with  the  added  ability  to 
prescribe  known  nodal  force  at  each  time  increment. 

The  dynamic  finite  element  code  with  the  updated  fracture  mechanic 
package  was  used  in  its  generation  mode  to  calculate  the  dynamic  fracture 
toughness,  Kj^.  Figure  1 shows  a typical  finite  element  breakdown  of  a M-CT 
specimen  used  in  this  study.  The  prescribed  crack  tip  motions  used  in  these 
series  of  generation  calculations  were  obtained  either  from  the  streaking 
photographs  of  two  Homalite-100  and  one  polycarbonate  M-CT  specimens  or  from 
the  discrete  crack  tip  recordings  from  the  sixteen  photographs  of  two  poly- 
carbonate M-CT  specimens.  Also  the  past  practice  [8,9]  of  using  dynamic 
elastic  modulus  in  stress  wave  propagation  was  discarded  and  the  static 
elastic  modulus  was  used  throughout  all  static  and  dynamic  numerical  analyses. 
This  procedure  follows  the  conclusion  of  Schirrer  [24]  who  concluded  that  the 
variations  in  elastic  modulus  did  not  affect  the  stress  distribution  appreciably 
but  did  change  the  strain  distribution  around  the  propagating  crack  tip. 

Static  mechanical  properties  of  Homalite-100  and  polycarbonate  specimens  were 
obtained  from  References  [8,9],  respectively. 

RESULTS 

The  dynamic  fracture  toughness,  Kjp,  during  crack  propagation  and  arrest 
in  one  Homalite-100  and  one  polycarbonate  M-CT  specimen  are  shown  in  Figures 
4 and  5,  respectively.  Kj^  results  in  these  figures  were  generated  numerically 
from  the  crack  tip  motion  obtained  from  streaking  photography.  Also  shown 
in  these  two  figures  are  the  corresponding  static  stress  intensity  factors 
obtained  by  static  finite  element  analysis.  The  static  Kj  results  differ 
with  those  reported  in  Reference  [1],  due  to  difference  in  modelling  the 
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applled  load  In  the  M-CT  specimen,  and  Is  In  agreement  with  the  corresponding 

i 

results  obtained  from  compliance  calibration  [10].  K^p  In  Homa11te-100  and 
polycarbonate  M-CT  specimens,  both  precipitously  drop  and  continue  to  remain 
at  nearly  constant  Kjp  thereafter.  The  higher  static  Kj  with  respect  to  the 
dynamic  Kjp  Is  an  Indication  that  much  of  the  released  energy  during  crack 
propagation  In  these  specimens  Is  dissipated  through  kinetic  energy  without 
being  returned  to  the  crack  tip  for  dissipation  through  fracture  energy. 

The  results  of  the  Homallte-100  M-CT  specimen  In  Figure  4 Is  In 
qualitative  agreement  with  the  dynamic  photoelastic  results  obtained  for  a 
slightly  larger  Homalite-100  M-CT  specimen  (of  12.7  mm  thickness)  In  Refer- 
ence [21].  The  low  crack  arrest  toughness,  0.25  (227  psl/frT)  «> 

was  approximately  equal  to  the  minimum  dynamic  fracture  toughness  in 
Figure  13  of  Reference  [8]  and  Is  55  percent  of  the  fracture  toughness, 

Kj£  f 0.42  MPa^T”  (380  psi  The  crack  arrest  stress  Intensity  factor 

of  Kj^  = 1.65  MP^/ir  (1500  ks1  /TrT)  for  the  polycarbonate  M-CT  specimen 
In  Figure  5 Is  nearly  one  half  of  the  pop-in  fracture  toughness  of 
Kjj.  = 3.4  MPa  /m"  (3.1  ks1  /in’). 

Figures  7 and  8 show  the  Kjq  versus  a relations  obtained  by  dynamic 
photoelasticity  for  two  polycarbonate  M-CT  specimens.  Also  shown  are  the 
Kjq  generated  numerically  using  the  measured  crack  velocities  shown  in 
Figure  9.  As  discussed  In  Reference  [22],  slight  Irregularities  In  crack 
velocity  variations  with  crack  extension  contributed  to  the  modest  differ- 
ences In  computed  and  experimentally  determined  dynamic  fracture  toughness. 

Considering  the  Idealized  elasto-dynamic  model  used  In  the  dynamic  finite 
element  analysis,  the  agreements  between  the  experimental  and  numerical 
results  are  good. 
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DISCUSSIONS 

The  differences  between  the  static  and  the  dynatni c stress  intensity 
factors,  as  shown  in  Figures  4,  6,  7 and  8 are  substantial  and  do  not 
exhibit  the  quasi-static  response  observed  in  the  AISI  1018  M-CT  specimen 
of  Reference  [1].  The  Homalite-100  M-CT  specimens  in  Reference  [7]  and 
in  this  paper,  the  polycarbonate  M-CT  specimens  of  this  paper  and  A533B 
M-CT  specimens  of  Reference  [6]  all  exhibit  a characteristic  decrease  in 
Kjq  followed  by  a relatively  stationary  for  almost  3/4  of  the  crack 
extension  prior  to  crack  arrest.  The  Kjp  in  the  mild  steel  specimen  of 
Reference  [1],  on  the  other  hand  decreased  nearly  monotonically  and 
closely  followed  the  corresponding  static  stress  intensity  factor.  This 
discrepancy  between  the  bulk  data  and  the  singular  data  of  Reference 
[1]  could  be  attributed  in  part  to  the  highly  localized  brittle  weld 
starter  crack  used  in  the  latter.  The  artificially  low  Kjg,  possibly  lower 
than  the  of  mild  steel,  of  brittle  weld  crack  starter  initiated  crack 
propagation  under  unrealistic  low  static  fracture  toughness  and  was  followed 
by  a quasi-static  crack  propagation.  Since  the  purpose  of  developing  this 
M-CT  specimen  is  for  dynamic  testing  [10],  it  can  be  concluded  that  the 
AISI-1018  M-CT  specimen  with  a brittle  weld  starter  crack  did  not  fulfill 
its  intended  use. 

Figure  10  shows  the  Kjp  versus  a relation  obtained  from  the  results 
of  Figures  5,  7 and  8.  Superposed  on  this  figure  is  the  averaged  Kjp  versus 
a relation  obtained  from  polycarbonate  dynamic  tear  (DT)  specimens  [9].  Not 
only  did  the  Kjq  versus  a relation  for  the  M-CT  specimen  shift  slightly  towards 
the  lower  a but  appears  to  reach  a lower  maximum  crack  velocity.  The  latter 
result  is  consistent  with  the  corresponding  results  for  Homalite-100  specimens 
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[5].  On  the  other  hand,  the  Kj^  versus  a relation,  which  was  obtained  from 
Figure  4 for  the  Homalite-100  M-CT  specimen,  coincided  with  that  for  Homalite- 
100  DT  specimens  [8].  The  slight  difference  in  Kj^  versus  a relations  in 
polycarbonate  fracture  specimens  is  in  agreement  with  the  results  of  Kalthoff 
et  al.  [25],  which  are  also  verified  numerically  by  Hodulak  et  al.  [26]. 

Figures  11  and  12  show  the  computed  energy  partition  in  a fracturing 
Homalite-100  and  polycarbonate  M-CT  specimens,  respectively.  The  substantial 
kinetic  energy  term  in  the  polycarbonate  M-CT  specimen  in  contrast  to  that 
in  the  Homalite-100  specimen  and  is  probably  due  to  the  blunt  starter  crack 
used  in  the  former.  The  computed  energy  for  these  specimens  in  Figures  11 
and  12  as  well  as  other  specimens  balance  to  within  5 percent  of  total  input 
energy  as  an  indication  of  the  numerical  accuracy  which  can  be  expected  in 
these  analyses. 

CONCLUSIONS 

1.  The  Kjjj  variations  with  crack  propagation  in  the  limited  number  of 
Homalite-100  and  polycarbonate  M-CT  specimens  analyzed  are  consistent 
with  similar  findings  by  others  [6-9]. 

2.  The  significant  difference  in  dynamic  responses  between  the  mild-steel 
M-CT  specimen  with  brittle  weld  starter  crack  and  the  Homalite-100  and 
polycarbonate  M-CT  specimens  requires  further  investigation. 

3.  The  shift  between  the  Kjq  versus  a relations  between  Homalite-100 
M-CT  specimens  and  DT  specimens  as  well  as  in  polycarbonate  specimens 
could  be  another  indication  of  the  geometry  and  size  dependence  of  the 
Kjp  versus  a relation. 
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FIGURE  3.  STREAKING  PHOTOGRAPH  OF  A PROPAGATING  CRACK 
TIP  IN  A POLYCARBONATE  M-CT  SPECIMEN  (THICK- 
NESS 6.4mm). 
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FIGURE  4.  CRACK  VELOCITY  AND  STRESS  INTENSITY  FACTORS 
OF  A FRACTURING  HOMALITE-lOO  M-CT  SPECIMEN. 
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FIGURE  5.  DYNAMIC  FRACTURE  TOUGHNESS  AND  CRACK  VELOCITY 
OF  A FRACTURING  POLYCARBONATE  M-CT  SPECIMEN. 
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FIGURE  7.  STRESS  INTENSITY  FACTORS  OF  A FRACTURING  POLY- 
CARBONATE M-CT  SPECIMEN,  J- 280678. 
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FIGURE  8.  STRESS  INTENSITY  FACTORS  OF  A FRACTURING  POLY- 
CARBONATE M-CT  SPECIMEN,  J-020878. 
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FIGURE  10.  DYNAMIC  FRACTURE  TOUGHNESS  VERSUS 
CRACK  VELOCITY  RELATION  FOR  POLY- 
CARBONATE SPECIMEN. 


FIGURE  II.  ENERGY  OF  A FRACTURING  HOMALITE  - 100 
M-CT  SPECIMEN,  SI00277. 


FIGURE  12.  ENERGY  OF  A FRACTURING  POLYCARBONATE  M-CT  SPECIMEN, 
J -020878. 
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